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ABSTRACT 

We report the detection of pulsations with ~ 124 s period in V2116 Oph, the optical 
counterpart of the low-mass X-ray binary GXl+4. The pulsations are sinusoidal with 
modulation amplitude of up to 4% in blue light and were observed in ten different 
observing sessions during 1996 April- August using a CCD photometer at the 1.6-m 
and 0.6-m telescopes of Laboratorio Nacional de Astrofisica, in Brazil. The pulsations 
were also observed with the UBVRI fast photometer. With only one exception the 
observed optical periods are consistent with those observed by the BATSE instrument 
on board the Compton Gamma Ray Observatory at the same epoch. There is a definite 
correlation between the observability of pulsations and the optical brightness of the 
system: V2116 Oph had R magnitude in the range 15.3—15.5 when the pulsed signal was 
detected, and R = 16.0 — 17.7 when no pulsations were present. The discovery makes 
GXH-4 only the third of ~ 35 accretion-powered X-ray pulsars to be firmly detected 
as a pulsating source in the optical. The presence of flickering and pulsations in V2116 
Oph adds strong evidence for an accretion disk scenario in this system. The absolute 
magnitude of the pulsed component on 1996 May 27 is estimated to be My ~ —1.5. The 
implied dimensions for the emitting region are l.li?0, 3.2i?0, and 7.0i?Q, for black- 
body spectral distributions with T = 10^ K, 2 x 10^ K, and 1 x 10^ K, respectively. 

Subject headings: stars: individual (V2116 Oph) — pulsars: individual (0X1+4) — 
stars: oscillations — stars: neutron — X-rays: stars 



Based on data collected at CNPq/Laboratorio Nacional de Astrofisica, Brazil 
Also Instituto Nacional de Pesquisas Espaciais, Brazil 



1. Introduction 

GX 1+4 is a low-mass X-ray binary (LMXB) 
well known for its unique characteristics such as 
extended high/low states (McCIintock &: Lcv- 
enthal 1989; Mony et al. 1991; Makishima et 
al. 1988), luminosity of up to ~ 10^^ erg/s in 
the high state (Rao et al. 1994; Nagase 1989), 
the hardest spectrum of all X-ray binaries (with 
kT ~ 35 — 55 keV in a thermal bremsstrahlung 
model) (Laurent et al. 1989; Staubert et al. 
1995; Predehl, Priedrich &; Staubert 1995) and 
the fastest recorded rate of change of period 
amongst the persistent X-ray pulsars {P/P ~ 
-2%/yr in the 1970s) (Nagase 1989). The optical 
counterpart of GX 1+4 was identified by Glass Sz 
Feast (1973) and confirmed recently by an im- 
proved X-ray position obtained with the ROSAT 
satellite (Predehl et al. 1995). Davidsen et al. 
(1977) showed that V2116 Oph has many char- 
acteristics of a symbiotic system, with the late- 
type star being classified as M6 III. Shahbaz et 
al. (1996) classify the late-type spectrum as M5 
III. The absorption spectrum of V2116 Oph red- 
ward of Ha is very similar to the spectrum of 
NSV 11776 (Cieshnski, Elizalde & Steiner 1994) 
which is classified as M4 III or later. Davidsen 
et al. (1977) point out a well defined excess of 
blue light in the spectrum, and that Ha, as well 
as the (blue) continuum, show variability. 

Recent low-time-resolution photometric stud- 
ies have shown interesting results like the ap- 
proximately simultaneous Ha and X-ray flaring 
activity (Manchanda et al. 1995; Greenhill et 
al. 1995) and the very low levels of Ha emis- 
sion (Sood et al. 1996) following the X-ray low- 
state in which GX 1+4 entered on 1996 August 
18 (Chakrabarty & Prince 1996). 

Despite all of the interesting features shown by 
the X-ray source, surprisingly little optical pho- 
tometry work has been done with good time res- 
olution. In the sole short-term variability study 
of this system performed so far (Krzeminski & 
Priedhorski 1978), a search for optical pulsations 



in Ha established upper limits of 1.7% and 0.7% 
to the pulsed fraction of sinusoidal and X-ray 
pulse shapes in the range 111 s < P < 143 s. This 
null result was interpreted as due to smearing and 
suppression effects due to the large dimension of 
the Ha emitting region, rua ~ 6 x 10^^ cm, as 
calculated by Davidsen et al. (1977). The es- 
cape time for Ha photons with an optical depth 
THa ~ 10^ from a cloud of that size would be 
(rHa/cjTHa ~ 2 X 10^ sec. 

Motivated by the discovery of flickering with 
time scales of ~minutes (Braga et al. 1993), and 
night-to-night variations on the mean R mag- 
nitude of up to 0.5 mag, we are conducting a 
concentrated observational effort trying to an- 
swer fundamental questions about this system: 
what is its orbital period? is this a wind-fed or 
accretion-disk-fed pulsar? What is the origin of 
the photometric variability in long and short time 
scales? 

2. Observations and Data Reduction 

The data discussed in this paper were ob- 
tained with the 1.6-m and 0.6-m telescopes of 
Laboratorio Nacional de Astroffsica, in Brazil. 
The CCD photometer we have built consists of 
a Wright Instr. thermoelectrically-cooled cam- 
era with a EEV CCD-02-06 back-illuminated 
chip operated in the frame-transfer mode (385 
X 289 pixels in the image section). The tim- 
ing for the instrument is provided by a Datum 
Inc. BC627AT Global Position System (GPS) 
receiver and timing board installed in the same 
IBM-PC clone used to read the CCD camera. 
Integration times down to 1 second are possible, 
especially with 2x2 on-chip binning. 

The plate scale produced by the various com- 
binations of telescopes, focal reducers and bin- 
ning factors is always ^ l'.'2/pixel. The field 
around V2116 Oph is relatively uncrowded, per- 
mitting safe aperture photometry to monitor its 
brightness relative to several comparison stars. 

The reduction of the data is presently done 
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off-line using a set of IRAF|^ scripts written by 
one of us (FJ). The reduction of a few thousand 
images, including the standard preparation pro- 
cedures (bias, dark and flat-field) can be done in 
the day following the observations. 

The data summarized in Table 1 resulted from 
20,465 images in the field of V2116 Oph. The 
CUSO4 filter used is a 4mm-thickness liquid solu- 
tion similar to what we have used for many years 
in the U passband of our UBVRI photometer 
(Jablonski et al. 1994). This non-standard pass- 
band maximizes photon count without sacrific- 
ing too much the capability of doing differential 
photometry. Column 3 of Table 1 shows the in- 
tegration time At, and column 4 shows the total 
duration of each run. On 1996 April 26 for in- 
stance, the data set consists of 3300 sequential 
images of 5 s integration time each. 

Each night, a few images in the R filter were 
obtained to provide a long-term light curve of 
the system. Since we have determined photoelec- 
trically the magnitude of star ^9 in the chart 
of Doxsey et al. (1977) {V = 13.39 ± 0.04, 
V - R = +0.77 ± 0.02, R - I = 0.79 ± 0.01) 
even nights of poorer photometric quality pro- 
vide good magnitude measurements for V2116 
Oph. This is the same comparison star used in 
the Ha photometry of Greenhill et al. (1995). 
Column 5 of Table 1 lists the mean R magnitude 
of V2116 Oph in each observational run. 

3. Analysis 

Figure 1 shows the light curve of V2116 Oph 
obtained with the CUSO4 filter on 1996 April 26, 
where the optical pulsations with 124 s were first 
detected (Jablonski et al. 1996). The compari- 
son stars observed simultaneously (two of them 
are shown in this figure) provide a very good es- 
timate of the noise in the variable star's light 

■^IRAF is distributed by National Optical Astronomy Ob- 
servatories, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under contract 
with the National Science Foundation. 



curve. 

Figure 2 shows the power spectrum of the light 
curve in Figure 1. The dotted line marks the 
expected level of noise estimated from the light 
curves of six comparison stars measured simulta- 
neously. The isolated feature close to 0.008 Hz 
has peak power ^ 200 times the power in the 
local continuum. A discrete periodogram with 
increased frequency resolution gives a better es- 
timate of the period of that feature: 124.17ib0.04 
s. There are two interesting additional pieces of 
information in Figure 2: the inclination of the 
power spectrum at low frequencies, a = 2.01 it 
0.05 in Power oc /~", and the period where this 
component intercepts the noise level, Pq = 57 it 2 
s. The flickering behavior of V2116 0ph is very 
similar to that shown by cataclysmic variables 
(Bruch 1992). An estimate of the amount of flick- 
ering in the light curves is given in the last col- 
umn of Table 1. It measures the excess power 
above the photon noise level, integrated in the 
frequency range 0.001 — 0.005 Hz. We follow van 
der Klis (1989) and express this quantity in terms 
of percentage r.m.s. variation in the light curve. 

The detection of stable optical pulsations with 
period close to what is expected from the X-rays 
history of GX 1-1-4 prompted us to observe the 
system as often as possible. Table 1 shows a sum- 
mary of all subsequent observations. In the case 
of a detection, the barycentric period and the 
semi-amplitude of the pulsed signal are shown in 
columns 6 and 7. In Figure 3 we plot the pe- 
riods listed in Table 1 together with the results 
from the continuous monitoring by BATSE. The 
very good quality of the X-rays daily period de- 
terminations results from the large modulation 
fraction and from the use of five-day segments in 
the analysis (Chakrabarty et al. 1997). 

4. Discussion 

We can summarize the results of this paper as 
follows: 

We have detected optical pulsations in V2116 
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Oph with periods very close to the expected from 
the X-rays history of GX 1+4 (Chakrabarty et 
al. 1997). This adds a third object to the smah 
group of accreting pulsars that are LMXB and 
present optical pulsations. The other cases are 
HZ Her=Hcr-Xl (Davidson et al. 1972; Mid- 
dleditch & Nelson 1976) and KZ TrA=4U1626-67 
(Ilovaisky, Motch & Chevalier 1978; Middleditch 
et al. 1981). 

We interpret the presence of flickering in the 
light curves as evidence of an accretion disk in the 
system. The power spectrum shown in Figure 2 
indicates that the variations due to flickering ap- 
pear on time scales as short as 57 s before they 
get immersed on the flat noise component due 
to photon statistics. If we associate the shortest 
time scale for flickering in our data, namely 57 s, 
with the radius of a Keplerian orbit in the accre- 
tion disk, we obtain a dimension r ~ 2.5 x 10^ 
cm, for a 1.4Mq neutron star. In other words, if 
we associate r with an upper limit on the magne- 
tospheric radius tm (Prank, King & Rainc 1992), 
then the fastness parameter uig in the theory of 
Ghosh & Lamb (1979a,b) would be « 0.5. This 
is consistent with a;^ < 1 needed for steady ac- 
cretion. The standard theory (Frank et al. 1992) 
also gives an estimate of the surface magnetic 
field of the neutron star, Bq 10^^ G, assum- 
ing a typical luminosity of 10^^ erg This 
is consistent with a previous estimate by White 
(1988). 

The absence of flickering when the system is 
faint (the total r.m.s. variation above photon 
noise in the light curve of June 19 is ^ 1%) in- 
dicates that the accretion may stop episodically. 
This is consistent with the observed low states 
in X-rays (Chakrabarty k, Prince 1996) contem- 
poraneous with very low levels of Ha emission 
(Sood ct al. 1996). If Roche lobe overflow is 
the preferred accretion mode in GXl-l-4, we are 
forced to conclude that the M companion does 
not fill its lobe all the time. This is expected 
if the M-giant presents large-scale mass motions 
it its outer layers, like red irregular and semi- 



regular variables usually do (Querci 1986). Be- 
sides providing episodes of low accretion, semi- 
regular variations in the accretion flux could also 
be a good explanation for the so called "300 
days" cycle in the torque history of the neutron 
star (Cutler, Dennis k, Dolan 1986; Chakrabarty 
et al. 1997). 

The origin of the optical pulsed emission can- 
not be determined from the data discussed in 
this paper, but taking what happens in HZ Her 
and KZ TrA as a paradigm, it could be a com- 
bination of reprocessing of X-rays in the accre- 
tion disk, in the atmosphere of late-type star, 
and in the stream between the latter and the 
disk (Middleditch k Nelson 1976; Middleditch 
et al. 1981). From the observed colors of the DC 
and pulsed light on 1996 May 27, we estimate 
the absolute magnitude of the pulsed component 
to be My ~ —1.5 for a distance of ~ 10 kpc 
and E{B — V) in the range 1.8 — 2.0 (Jablon- 
ski k Pereira 1997). If the pulsed light follows 
a black-body spectral distribution with temper- 
ature Tbb = 10^ K, the size of the emitting re- 
gion would be ~ l.li?Q. For lower temperatures, 
2 X 10^ K and 1 x 10^ K, we obtain dimensions of 
3.2i?Q and 7.0Rq, respectively. The correspon- 
dent areas are not prohibitive either as fractions 
of the surface area of the companion star or as 
fractions of an accretion disk area. The accretion 
stream is ruled out as a likely sole source of the 
pulsations because its projected area is too small. 
If the pulsations show a dependence with orbital 
phase, like in HZ Her and KZ TrA, the long term 
monitoring of the optical pulsations may be use- 
ful to answer a long-lasting (and very basic) ques- 
tion about this system: what is the orbital period 
of GX 1-1-4 ? So far, we only know that in order 
to fill the corresponding Roche lobe, a M4-6 III 
star with 80-100 Rq would need Porb ^ 100 days 
in a binary with mass ratio g ~ 1. 

The BATSE measurements (Chakrabarty et 
al. 1997) show that it is very difficult to measure 
the orbital period of the system from the Doppler 
effect on the X-ray pulses alone, by the follow- 
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ing reasons: the orbital period is long, proba- 
bly ^ 100 d, the orbital inclination may be low 
(though this is less certain), and it is very dif- 
ficult to separate the effects of the spin-down 
(and its fluctuations) from the Dopplcr effect it- 
self on a long-period orbit. Our best determi- 
nation of pulse period, namely 124.17±0.04 s, 
on April 26, is significantly different from the 
BATSE measurement, 123.9453 ±0.0014 s, based 
on a linear interpolation between the daily val- 
ues of Chakrabarty et al. (1997). If caused by 
Doppler effect, the implied difference in velocity 
between the X-ray and optical emitting regions 
would be ~ 540(sin)^^z km/s, too high for a bi- 
nary with the characteristics discussed earlier. 

By analogy with HZ Her and KZTrA we ex- 
pect to see sidebands of Popt or Px according to 
-^opt = ~ -^orb- explain the ob- 

served difference between Popt and Px in terms 
of sidebands, then the derived Porb is too small, 
~ 0.8 days, again inconsistent with our previous 
knowledge about this system. A good observa- 
tional coverage will settle this point, in particular 
if we are able to track the phase of the pulsations 
from simultaneous measurements in the optical 
and X-ray bands. 

The relatively long period and large ampli- 
tude of the optical pulsations in V2116 Oph open 
interesting perspectives for future work. Likely 
projects include the detailed history of the pul- 
sations on a night-to-night basis and its correla- 
tion with the X-ray timings, a spectrophotomet- 
ric study of the optical and NIR spectrum of the 
pulsations and the correlation of the visibility of 
the pulsations with orbital phase. 

MP acknowledges the Brazilian agencies CNPq 
and CAPES for graduate studies support; FJ and 
JB acknowledge CNPq for research fellowships. 
We thank Dr. Dccpto Chakrabarty for calling 
our attention to the incorrect error bar for the 
124 s period in Jablonski et al. (1996). The cor- 
rect value is shown in Table 1 of this paper. 
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Fig. 1. — Differential photometry CCD light 
curve of V2116 Oph on 1996 April 26, with 5 
s time resolution. The differential light curves 
of two nearby stars are shown for comparison. 
The brightest star is object 7^ 10 in the chart 
of Doxsey et. al (1977). The inset shows a sec- 
tion where individual pulses can be seen. The 
arrows pointing to subsequent maxima were pro- 
duced by constant-period (124.17 s) prediction. 
The small bar in the upper right corner of the 
inset shows the estimated level of noise in that 
section of the light curve. 

Fig. 2. — The power spectrum of the light curve 
in Figure 1. The horizontal dashed line shows 
the estimated level of noise derived from the light 
curves of six comparison stars measured simulta- 
neously. The inset shows the light curve folded 
on the 124.17 s period. The phase interval is 
repeated twice for better visualization and the 
vertical scale is relative amplitude. 

Fig. 3. — The period of the optical pulsa- 
tions listed in Table 1 (filled circles) together 
with the BATSE daily measurements (crosses; 
Chakrabarty et al. 1997). 
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Table 1 

Fast Photometry of V2116 Oph 



Date 


Filter 


At 


Duration 


R 


Barycentric period 


Semi-amplitude 


Flickering^ 


(1996) 




(s) 


(hrs) 


(mag) 


(s) 


(%) 


(%) 


Apr 26^ 


CuS04 


5 


4.6 


15.48 


124.17±0.04 


1.3±0.1 


2.3 


May 21^ 


CuS04 


15 


4.2 


15.26 


124.1±0.3 


0.6±0.3 


2.5 


May 25^= 


CuS04 


20 


0.4 


15.29 


124.6±1.6 


1.8±0.5 


0.7 


May 26^= 


CuS04 


10 


0.9 


15.29 


124.2±0.9 


1.2±0.4 


1.9 


May 2T 


CuS04 


20 


0.9 


15.28 


124.0±0.3 


4.3±0.5 


1.7 


May 27^= 


R 


20 


0.8 


15.28 


124.1±0.4 


2.6±0.4 


2.0 


May 27^'^ 


Clear 


15^ 


0.7 


15.28 


124.5±0.5 


2.5±0.4 


1.4 


May 27^"^ 


U 


15^ 


0.7 


15.28 




<12 


8.6 


May 27^ 


B 


15^ 


0.7 


15.28 




<4.4 


4.8 


May 27^"^ 


V 


15^ 


0.7 


15.28 


124.1±1.0 


5.3±1.6 


4.6 


May 27^^^ 


R 


15^ 


0.7 


15.28 


124.6±0.6 


3.3±0.6 


1.7 


May 27^^^ 


I 


15^ 


0.7 


15.28 


124.1±0.6 


2.2±0.4 


1.3 


TV It ^ OOC 

May 28 


UubL)4 


20 


1.5 


15.37 


123.9±0.5 


1.7±U.D 


2.5 


May 28^= 


R 


20 


1.7 


15.37 


124.3±0.5 


0.6±0.2 


1.1 


Jun 19^ 


Clear 


5 


7.0 


17.72 




<0.05 


0.1 


Jun 20^" 


Clear 


5 


6.6 


17.68 




<0.07 


0.2 


Jul 13^ 


R 


25 


3.7 


17.01 


124.1±0.2 


0.9±0.3 


1.2 


Aug 04^ 


R 


20 


3.0 


15.67 


124.2±0.3 


0.4±0.2 


0.5 


Aug 05^ 


R 


20 


4.4 


15.45 


124.0±0.4 


0.2±0.2 


1.3 


Aug 06^^ 


CuS04 


20 


3.3 


15.28 


124.7±0.3 


0.4±0.1 


0.6 


Aug 07*^ 


R 


20 


1.7 


15.29 




<0.1 


0.5 


Aug 18^ 


Clear 


10 


1.4 


15.99 




<0.15 


0.2 


Aug 19^ 


R 


20 


3.0 


16.29 




<0.33 


0.5 


Aug 20^^ 


R 


25 


4.2 


16.09 




<0.17 


0.3 


Aug 22"^ 


R 


25 


5.7 


16.94 




<0.16 


0.9 



^1.6-m telescope 

^0.6-m Zeiss telescope 

^0.6-m Boiler & Chivens telescope 

'^FOTRAP photometer 

®Time resolution instead of integration time 

^Integrated in the range 0.001 — 0.005 Hz 
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